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I Abstract 

Oh 

^ We construct an M-theory background dual to the metastable state recently discussed 

by Klebanov and Pufu, which corresponds to placing a stack of anti-M2 branes at the 
tip of a warped Stenzel space. With this purpose we analytically solve for the linearized 

>- 

non-supersymmetric deformations around the warped Stenzel space, preserving the 50(5) 
symmetries of the supersymmetric background, and which interpolate between the IR and 
CN UV region. We identify the supergravity solution which corresponds to a stack of N back- 

reacting anti-M2 branes by fixing all the 12 integration constants in terms of N. While in 
the UV this solution has the desired features to describe the conjectured metastable state of 
the dual (2+l)-dimensional theory, in the IR it suffers from a singularity in the four-form 
flux, which we describe in some details. 
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1 Introduction 

The work of Intriligator, Seiberg and Shih [T] has drawn attention to mechanisms of metastable 
supersymmetry breaking in quantum field theories. Since the constructions of such states in- 
volve strongly coupled regimes, it is natural to address the study of this phenomenon in stringy 
realizations of the supersymmetric theories and indeed there exist many corners of the string 
theory web where these constructions have been proposed. In Type llA string theory one can 
try to build D-brane models in order to reproduce metastability in the field theories engineered 
on the brane world volumes; however, it turns out that the probe brane picture (i.e. Qs = 0) is 
too naive and once the backreaction of the branes is taken into account these systems can fail to 
reproduce metastable states of the supersymmetric theories [2]. 
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Another approach is to consider brane reahsations which extend the AdS/CFT correspon- 
dance to non-conformal or less-supersymmetric theories and try to construct metastable states 
in this context. One way to achieve this is to start by a configuration of branes placed at 
some Calabi-Yau singularity and consider the supergravity solution obtained after smoothing 
the singular point. The most well known example in Type IIB string theory is the Klebanov- 
Strassler (KS) background [3]. In this setting the first evidence of a metastable state in the 
SU {M{k + 1) — p) X SU {Mk — p) theory was given in [H |5], with a construction that involves 
placing a stack of anti-D3 branes in the KS geometry. These branes are attracted to the tip of 
the deformed conifold and polarize into an NS5-brane due to the Myers effect [6]. In a probe 
analysis, it was shown [4j that this state is metastable and long-lived if the number p of anti-D3 
branes is small (~ 8%) compared to the M units of R-R 3-form fiux of the unperturbed super- 
symmetric background. It is an important question whether this picture remains valid once the 
backreaction of the anti-D3 branes on the KS geometry is taken into account. Constructing the 
full backreacted supergravity solution is a difficult task, but one can make progress if some sim- 
plification is adopted, namely i) smearing the sources on the of the deformed conifold and ii) 
working in perturbation theory around the supersymmetric background, at first order in the pa- 
rameter p/M. This simplified (but still difficult enough) problem was investigated in (TIE], where 
the full linearized backreaction of the anti-D3 branes on the KS geometry has been constructecQ 

It is of obvious interest to address the same question in different contexts where metastable 
states are conjectured to exist, by string theory arguments, in supergravity backgrounds dual to 
strongly coupled field theories. 

In this paper we study the case of an AdS^^/CFT^ correspondence, which involves an A/" = 2 
supersymmetric (2+l)-dimensional theory, whose supergravity dual is AdS^^ x V^^2i where V^^2 
is the 7-dimensional Sasaki-Einstein space V^^2 = S0{5)/ S0{3). Recently, a gravity dual for a 
long-lived metastable state has been proposed in [10] based on the probe analysis, by placing a 
stack of anti-M2 branes at the tip of the warped M-theory background with transverse Stenzel 
space [llj. Here, the analogue of the KS solution is the supersymmetric solution of Cvetic, 
Gibbons, Lii and Pope (CGLP) [12]; indeed, the 8-dimensional Stenzel space is a part of a family 
of Ricci-fiat solutions parametrized by the dimension n, which include the deformed conifold for 
n = 6. The mechanism for which the false vacuum decays is similar to the KPV process [4j: the 
anti-M2 branes fall in the warped throat and at the tip they polarize into M5-branes wrapping 
an C 5*^. The probe analysis of (TU] shows that this state is metastable if p/M < 0.054, where 
M is the number of units of the 4-form fiux of the CGLP background. 

The effects of the backreaction of the anti-M2 branes on the transverse geometry have been 
studied in ^13j, where the linearized equations that govern the first-order backreacted solution 

-'^This problem was also adrcsscd in [Sj. 
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have been solved implicitly in terms of integrals by using the first-order formalism introduced 
by Borokhov and Gubser [Tl] and the full solution was presented separately in the small and 
large radius limit. The main purpose of that work was to study the IR behavior of the perturbed 
solution, and the conclusion of this analysis was similar to the anti-D3 case, namely that the 
conjectured solution dual to the metastable state exhibits certain singularities which in the anti- 
M2 case lead to a divergent action in the IR[^ To decide whether this singularity is admissible or 
not is a difficult task, and the answer is clearly beyond the linearized approximation. One way to 
proceed is to connect the IR and the UV region and to see if the conjectured solution eventually 
develops problems in the ultraviolet. In this perspective, it is clearly interesting to perform such 
an analysis in the anti-M2 brane configuration, which in the IR can be thought as the M-theory 
generalization of the Type IIB KS solution, but has a rather different behavior in the UV. For 
example, in the M-theory background there is not a logarithmic running of the charge, which 
is an important feature of the KS background, and was crucial in the analysis of the backreac- 
tion performed in [S|. In this paper we perform the analysis outlined above and by extending 
the results of [I3j we present the full analytic solution of the linearized supergravity equations 
which describe the most general non-supersymmetric deformation of the warped Stenzel space 
compatible with the symmetries of the CGLP background. With this result we are able to study 
the effects of IR boundary conditions on the ultraviolet behavior of the supergravity modes and 
we identify the unique solution which has the desired features to describe anti-M2 branes in the 
CGLP background (leaving open the issue of the singularity discussed above). 

This paper is organized as follows. In Section 2 we review the computational formalism 
and we solve analytically the system of first-order differential equations governing perturbations 
around the CGLP supersymmetric solution. Our full solution, which is shown in Appendix |A| 
contains few single integrals that cannot be explicitly performed, but they can easily be handled 
with numerical integration. In Section 3 we show the expansions of our solution in the IR and 
in the UV region in terms of a set of twelve integration constants denoted (X^,, Ya). In Section |4] 
we discuss the various charges in the Stenzel background and we identify the perturbation due 
to the presence of M2 branes. In Section |5] we impose the boundary conditions that arise from 
placing a stack of anti-M2 branes at the tip (r = 0) of the geometry and we discuss the problems 
associated to an infrared singularity in the fluxes. We then summarize the asymptotic behavior 
of the anti-M2 solution, which is expressed in terms of the number of anti-M2 branes. As a 
check of our boundary conditions, we compute the force exerted on a probe M2 brane and we 
show that it agrees with the one derived from the brane /ant ibrane potential (which we review 
in Appendix [B]) . We end with a discussion in Section 6. 

^See [TB] for a similar analysis in a Type IIA context. 
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2 Linearized equations and their solutions 

The linearized equations governing the deformations around the warped Stenzel space have been 
derived in [13] by using the Borokhov-Gubser [H] first-order formahsm. We use the ansatz for 
the S'0(5)-invariant supergravity solution of rrnisi 



where cxj, di {i = 1,2,3) and u are the 1-forms in the coset S0{5)/ S0{3) and /i = 0, 1,2. The 
four-form field strength G4 is given by 

G4 = dK{r) A dx^ A dx^ Adx"^ + mF4, (2.2) 

where the internal flux F4 is parametrized by 

F4 = f'{r)dr A ai A 5-2 A CTg + h'{r) e'^^dr A A A at (2.3) 



+ -{Ah{r) - f{r))e^^^u AdiAajAak-G h{r) u Aai A02 Aa^. 



1 

2' 

and the function K{r) is fixed in terms of the other functions by the equation of motion 



rf^G^ = AG4. (2.4) 

The method introduced in [Tl] relies on the existence of a superpotential W defined such that 
its square gives the potential, namely 

1 , dW dW 

We consider an expansion for the fields 0" (a = 1, ra) around the supersymmetric background 

r = r, + rAx) + o{x^), (2.6) 

where X represents the set of perturbation parameters, is linear in them, and (p^ are the 
functions in the CGLP solution, written explicitly in (2.15). We will denote the set of functions 
0", a = 1, 6 in the following order 

r = {a,[5,^,zj,h). (2.7) 

The first order formalism gives a set of 2n linear first-order differential equations for the pertur- 
bations (pl and their conjugates 

^ + 6M^(0o) = O, (2.8) 
'^'^^ M^(0o)0? = G"^''a, (2.9) 



dr 

■^A more general solution which includes this 50(5)-invariant ansatz has been constructed in [15j 
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where 



= G, 



dr 



M\,= 



1 d 



G 



be 



dW 



(2.10) 



The equations (2.9 ) are the definitions of the ^a,while the n equations (2.8 ) form a closed set and 
imply the equations of motion [H] . The functions should additionally satisfy the zero-energy 
condition 



dr 



0. 



The field-space metric in (2.10) is 

1 



G 



ah ' 



-a-3(/3+z) 



and the superpotential is given by 

W{<j)) = -3e2"+2/5(e2" + e^^ + e^^) - 6m^e-'' \h{f - 2h) 
The background fields satisfy the fiow equatioiiQ 



1 1 
54 



dr 



2 04 



and they are given by the CGLP solution [12 



,2ao ^ l(2 + cosh2r)^/^coshi 
3 



.2/30 
,270 

/o 
ho 



- (2 + cosh 2r) ^^"^ sinh r tanh i 
3 

(2 + cosh2r)"^/^ coshV , 



1 — 3 cosh r 

3^/^ cosh^ r ' 
1 

2 3^/^ cosh r ' 
- log(m^if(r)) , 



where the warp factor H is defined by the following integral: 
3 sech tanh u 



H(r] 



(2 + cosh 2m) 3/4 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



= ^ ( JJ_ 1^)3/2 + 5v/2 F ( arcsin(y-i) I - l) , (2.16) 



^Note that the equations for a = 1,2,3 are equivalent to the the Ricci-flat Kahler condition for the metric [17] . 
while the ones for a = 5, 6 give the self-duality of the internal form F4. 
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where 

y = (2 + cosh(2r))^/^ (2.17) 
and F is the incomplete elhptic integral of the first kind 

F{<p\q) = [\l - q sm\9))-'/^d9 . (2.18) 
Jo 

As shown in [13] , it is useful to solve for the following linear combinations of the fields and (pa 

L = (6 + 6 + ^3, ei - 6 + 3 6, 6 + 6 - 36, ^4, ^5, ee) , (2.19) 

<f>a = (01 - 4>2, 01 + 02-2 03, 03, 04, 05, 06) . (2.20) 

The first-order systems of coupled differential equations for the fields and 0a are: 

el 



and 



= ie-3(^o+"o+'^o)^2 ^^^^^ _ 2;,^) _ 1) 1^ ^ 
y 




(2.21) 


= ^e-=^(^o+°°+'^«)m2(54/io(/o - 2ho) - l) ^4 , 
y 




(2.22) 






(2.23) 


= 6e-^"°+^^«f5 - 2e°»-*f6 - 2m2(/o - Aho)e 


-3(20+00+^0) f 

S4 5 


(2.24) 


2 
9 


-2/io)/io-l)l4, 


(2.25) 


= 2cosh(«o - /3o)6 - ^e°°-'^«ei + ^e-°°-'^«( 


e2"o _ 2e2To)f3 


(2.26) 



36/ioe-=^"°+=''^°e5 + e"«-*(/o - 4/10)^6 



-2 cosh(ao - /3o)0i + ^e-3K+/3o)(_3|^ + + 3^3) ^ (2.27) 
_4g-ao-/3o+27o02 - 6 sinh(ao - /3o)0i + ^e-^("«+^°)(-3|i + 7^3) , (2.28) 



3 

2^ ■ 12" 

2 



.3 = 3 sinh(ao - /3o)0i + ^e-^o-^^+^^o^^ + _e"3K+/3o)(|^ _ 3^3) ^ (2.29) 



()5 = ^e-3°»+=^''« (-3m'06 + 9m'/io0i + e^^^^s) , (2.30) 
1 



e-o-/^o(_3^2^^ ^ ^2mV6 - 3m^(/o - 4/io)0i + e^^^^g) , (2.31) 



^1 = le-=^(^°+"«+^°) (2e3^«e4 + m\{l + 54/io(2/io - /o))02 + 203 + 04 (2.32) 
+ 18(/io(-3(/o - 2/lo)(203 + 04) + 05) + (/o - 4/lo)06))) • 



2.1 Solutions for 

The solution for the modes has been derived in [T3]. We first note that the equation for ^4 



can be easily integrated by using the fiow equations (2.14); for a = 4 this reads: 



This shows that ^4 is proportional to the warp factor: 

^4 = m^H{r)X^. 

In terms of the radial variable r the solutions for the remaining modes are 



^4 = m''X,Hir) , 
^5 = — eschar sechr ( ^ 



\ 2 
3v^ 



m^X^ 



csch 



(cosh2M + 2)3/4 

cosh 1u csch sech 
(cosh2M + 2)3/4 

3v^ r csch'^u 



, 2 / 3 V3 2 X. /■ 

— cosh r coth r ^ 2~"^ / 

^6 = (3 cosh 2r + 1) csch sech ^ ^ 771^X4 y 

(3 



-771^X4 



+ Xg j , 
(iw + X5 

du + Xe 



(cosh2M + 2)3/4 
cosh 2u csch sech 



(cosh 2m + 2)3/4 



-, A / , N / 9 ; csch sech , 
^3 = -6 sinh^ r (cosh 2r + 2) 771^X4 / — — + X; 



^2 = sinh r cosh r 



(cosh 2m + 2)7/4 

3 4 

X2 H — Xi coth r + 9X3 sinh^ r cosh r H csch r sech 

2 3v 3 



4 x^ / , o 1 \ /^oTT/ \ 1 2tanhrsech r 

H ^Xg (coth r — 3 tanh r j + 777 X4 3-n ( r ) coth r — -— r- 

3v3 ^ (cosh 2r + 2)3/4 



+ 2 csch r sech r 



(cosh 2m + 2)3/4 



(2.33) 



(2.34) 



(2.35) 



eschew o -, r csch sech '^■u 

raM + 9 smh rcoshr / ^-z ^,„,, du 



(cosh 2m + 2)7/4 



r' cosh 2m csch ^m sech ^m 

+ 4(cosh2r — 2) csch2r / ; ; -— ; du 

^ ' ' (cosh2M + 2)3/4 



The zero energy condition (2.11) reads 



Xo = 0. 



(2.36) 
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By using the change of variables (2.17) it is easy to show (see Appendix [A]) that the solution can 
be expressed in terms of only two integrals, the warp factor H{r) and the Green's function [18] 



G(r) 



Sv^csch^M ^ a/3 
2 (cosh 2m + 2)3/4 " ~ ^ 



1/2 



My' - 1) 

\/2(9-3?/4) 



3 2-^/2f( arcsin(y-^)|l 



- V2U(^ - ^3; -arcsin(?/~^)|l) - ^2 n(^y3; - arcsin(?/-^)|l) 



(2.37) 



where we use the standard definition for the incomplete elliptic integral of the third kind 

U{n;(j)\m) = {1 - nsm'^{e))-\l - msm'^{e))-^^^de . 
Jo 



(2.38) 



2.2 Solutions for 0^ 

We now present the solution for the (pa modes. Here we show the result in a compact form in 
terms of the variable r and we relegate to Appendix |A]the involved analytic expressions which are 
obtained by explicitly performing the integrations. The first-order perturbations to the metric 
modes and fluxes are 



^ 9 cothucschw ~ Yi 
4i23 du+ - 



sinh2r J 2 (2 + cosh 2m) 3/4 
9 csch r sinh u 



4 (2 + cosh 2r) J (2 + cosh 2m)3/4 
3 7 csch 



sinh 2r 

(15 + 3 cosh 2u) - 12 ^2 - (23 + 7 cosh 2u) fa 



(2.39) 

du 



+ 



2 + cosh 2r 2 + cosh 2r 
9 r csch '^u 



Yo 



32 J (2 + cosh2M)3/4 



^1 + 3 cosh 2u ^123 + 3 ^3 



du —"^ cosh 2r 01 — ^ 02 
8 8 



65 = sinh r tanh r A5 + - cosh r (5 — cosh 2r) Ag , 

1 1 o 

66 = ~ ^ (3 + cosh 2r) sinh r tanh r A5 + - cosh r Ag 
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where we defined 



A. 



+ 



24 



12 



12 sinh u ^^5 ~ (5 — cosh 2u) coth^ m csch u 

(2 — cosh2'u)csch'^M ~ , 2 — cosh2r ~ 
^123 du + 



Hiu] du 



(2.40) 



(2 + cosh 2m)3/4 



6^3 



1 ; 



sechwtanh u 



3 (3 + cosh 2u) ,^5 + cosh^ u 



Hiu) du 



(2 + cosh 2u) ^/^sech "^u tanh -u ^123 du 



(2 + cosh 2r) tanh^ r ~ 

* 



and we dubbed ^123 the following combination of 

^123 = 36-46-36. 

The last mode we solve for is the perturbation to the warp factor 04. Its integral expression is 
1 /"^ 6m^ csch iff ~ . 1 3m^ sech^utanhu (02 + 203) 



6m^csch^Mif(u) ~ 
m^E{r) J (2 + cosh2M)3/4 " + 



rn^Hir) 



(2 + cosh2M)3/4 



du 



3v^ csch u sech u (csch 05 + 2 sech ^ 



du + 



(2.41) 



mm{r) J (2 + cosh2M)3/4 ' m^H{r) 

We now briefly explain the procedure we followed in order to obtain this solution. We firstly 



solve the system (2.27)-(2.32) using the Lagrange method of variation of parameters. While in 



principle this produce a solution with an increasing number of nested integrations, we found that 



successive integrations by parts reduce the outcome of this method to the compact form (2.39). 



We note that since the solution for the 6 modes is analytic in the variable the aforementioned 
solution for the 0a modes contain at most single integrals of the from 



/(m)L(u) du, 



(2.42) 



where L(?/) is a combination of incomplete elliptic integrals and f{y) is a polynomial function of 
the variable y. In this form the expressions for the modes 0^ can be easily evaluated numerically, 
and thus provide a full interpolating solution which connects the IR and the UV region. 

The space of solutions we solved for is parametrized by twelve integration constants X^, Ya, of 



which only ten are physical since X2 can be eliminated through the zero energy condition (2.36) 



and Ys corresponds to a rescaling of the three-dimensional coordinates. In Appendix |X] we 
show the full solution obtained after replacing the analytic expressions for the modes 6 ^md by 
recasting some of the integrations in terms of incomplete elliptic integrals. We were not able to 
further simplify the resulting solution, but we stress that the crucial improvement that permits to 
easily handle numerical evaluation is the absence of nested integration (as opposed for example 
to what happens for the anti-D3 case [19]). 
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3 Asymptotic behavior 



In order to impose the desired boundary conditions we need to calculate the behavior of the 
solution presented in the previous section in the small and large r limits. For that we need the 
expansions of the elliptic integrals that enter in the expressions for the (pa modes. In the IR the 
first terms of the relevant functions are 

F(^csi„fe-)| - l) = F„ - ^ + + Oir^) , (3.1) 

2 4 

n( - V3; - a.csi„(,-)l - l) = + ^ - + 0(r«) . (3,2) 

n(V3; - ^cMv-')\ - l) - K. + - + + 0(,.«) . (3,3) 

where in order to keep notation intelligible, we used the following abbreviations: 

Fq = F(|arcsin (^7^) I - l) ~ 0.7896 , Ki = u(^ - VS;- arcsin (^7^) I - l) ~ -0.6142 
and K2 ~ —0.9102. We also encounter the constant 

2^ ^ 



K{-l)=2\l-Ti-] ^1.3110^ 



where K is the complete elliptic integral of the first kincQ Finally, we need the expansion of the 



warp factor (2.16) which is 



1 T 4 

H{r) = Ho-- 31/^2 + ^ + 0(r6) , Ho = -4 3^/^ + Sv^Fq ^ 0.3187 . (3.4) 

With these expansions we can easily find the IR behavior of the solution. In order to match with 
the UV behavior we only need to perform a numerical integration to find the expansions of the 
integrals that appear as the coefficients of X4 in the solution shown in Appendix |Aj 

3.1 Numerical matching 

We now briefly describe the numerical method used to relate the UV and IR expansions of 



the integrals that appear in the solution for the (pa modes. They are of the form (2.42), thus 



by using (3.1)-(3.3) we easily get the IR expansions for the integrands. By performing an 
indefinite integration we therefore get the desired expansions up to an integration constant which 
is generically different in the IR and in the UV. Since these integrals are divergent in the small 



^Defined as K{q) ^ F(f |g). 
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6 - 
4 - 
2 - 




Figure 1: The solution for the mode 06, for X2 = 0, Xi = X3 = X5 = = 1, X4 = 10, = 1, 
m = 1 (underlying blue solid line). The red and orange dashed curves correspond to the IR and 
UV expansions (respectively up to 20 and 15 terms). 



r limit but vanish at r = 00, we chose to do the definite integration in the range [r, 00]; in this 
way the UV integration constant is zero and we only need to match in the IR. This can be done 
up to an arbitrary precision p by fixing an tq smaller than the radius of convergence of the IR 
series, evaluating the IR expansion S of the indefinite integral at tq up to the appropriate order 
n and then fixing a constant k such that 



Sn{ro) + k-j f{u)L{u)du 

J rn 



< w 



(3.5) 



We kept a precision of p ~ 10, which we found enough for our purposes. As a check, we can 
verify that the expansions obtained in the aforementioned way approximates well the numerical 
solution for small and large r, as shown in Figure [T] for one of the perturbation modes. 



3.2 Infrared expansions 

We now show the IR expansions of the modes 0a, focusing on the singular behavior which is 
needed in order to impose boundary conditions in Section |5} These expansions appeared already 
in [13] and apart from making sure our results are fine in the present paper we relate Y^^ to 
y^^, a crucial step in order to try and write the backreacted state at linearized order for all 
radii. Here the integration constants Xa and Ya are those appearing in the analytic solution 
shown in section 



2.1 



and 2.2 and we defined the Ya as 



Ya = Ya+m'X^k^ 



(3.6) 



11 



k^^ = 7.45479 , k^^ = 0.301287 , k^^ = 0.112188 , k^^ = 0.576358 , k^^ = -0.00504419 . 

The constants fc^^ are obtained with the numerical procedure outhned in the previous subsection. 
We also impose the zero energy condition (2.36) and so in what follows we set X2 = 0. With 
these remarks and notations in mind, we now provide the IR expansions for the first-order 
perturbation modes 



1 



433/4^2 
1 



+ 



12r 



27Xi - 16V3{X5 + Xq) -30m'^ HoX^ 
QYi + y2 3^/^(-45y3Xi - 162v^X3 + 4OX5 + 112 X6)fs:(-1) 



(3.7) 



+ Y^-^ (l89 Xi - 80v^ (X5 + Xe) + 6 Xi{83 Ho - 3^/' 33)) + 0(? 



^ + -^Ho (63X1 + 432X3 - 7v^(5X5 + llXe)) + 



105 



35 



72 



33/4 Xfi 



+ 



^ 3V^ Xi - ^ 3 X3 + 2 3^/^ Xs + ^ 3^/^ Xe 
20 35 ^ ^ 5 ^ 



(3.^ 



„, 16^ 768X3 3 31/4 m2X4 16 (5X5 + 11X5) X 4F2 

+ iln — — Ai — ;= 

^ 5 ' 35 2 15^3 J 9 



1 

12^ 



6F1 - V2 3^/^(45v^Xi + 162^3 X3 -40X5 - 112X6)ir(-l) 



+ 



137|V4Xi ^ 12792 3V^X3 + 5^3 ^2 _ 1 (225 3^/^ X, + 9081 3^/^ X, - 205f.) 
25 175 675 



V 75 



2624 23 3^/^ m2 X4 164 (5 X5 + 1 1 Xg 

X3 — — — 

175 2 



+ 



Ho 
70 

2733/^X6 Y2 



- 63 Xi - 432 X3 + 7^3(5 X5 + 11 Xg)) - 



225v^ ^ ^ ^^^^ 
189 3^/^X1 1053 



40 



35 



31/^X3 



(3.9) 



+ 



2Ho 
105 



(63 Xi + 432 X3 - 7^3(5 X5 + 11 Xg 



^ , _ 1167 _ 67 _ 671 X.- 1033 

160 100 16 12031/4 



41 
360 



+ ^ ( - 10836 Xi - 69444 X3 + 6300 3^/^ X4 + 7^3 (725 X5 + 1847 Xg) J 



+ 



630^2 
4^2 



{21 Xi - 4y3(X5 + Xg)) + 30^2 3^/^ X4 Fq log ; 
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+ ^ (21 3^^ Xi - 210 VSm" - 4 3^/^ (X5 + Xg)) logr + e>(r) , 
05 = ^ - 12 3^/^ - 81 32/^ X3 + 18 3^/^ Xq + 10 Fe) (3. 10) 

- ^ (^21^3 Xi + 4 (27v^X3 + 10 X5 + 4 Xg)) + 0(7 

\/3//o"l'X4-2(X5 + X6) 



,^3 A 



I -^0 

12 



^ / 7v/3Xi 9^3X3 7m^X4 WX^ + X^ 
°V 80 20 2431/4 90 



H^m^X, 433=^/^X1 8133/^X3 , 25^5 + 106X6 , % , ^ 

^^V! ^ ^ 30 33/4 + y + ' (3-11) 

Sv4fz2 2^ 1 / 207v^Xi 2403v/3X3 54 ,3/4 



33 3i/4;^^_ 333 g,/4^^ 5X5-X, 



ffo V 40 

+ C(r°). (3.12) 



40 70 531/4 



3.3 Ultraviolet expansions 

Here we show the leading terms in the UV expansions of the perturbation modes (pa 
01 = -^^e-'-/' + 2e-^' Yi - e'^'/^ 2^1^(21 X^ + 81 X3 - IGv^Xg) (3.13) 
+ ^^e-^'-/^ (3267X3 - 1024v^X6) + 0(e-«'-) , 

'/•3-1^3-^ + ^^e -^^^^e +-e + 0(e ), (3.15) 

05 = ^IZ^e^'- - ^i^^^e'- + - 8^3 + 117^5 + r!Y^ (3.16) 

- 38 2^/4 V^e-^'-/^ X3 + ^e-^' ^168^3^1 - 9 (111 + Fe)) 

+ -|-2^/^e-^"/2(-3348\/3Xi + 1323\/3X3 + 6160X5) + Oie-^') , 
195 

06 = -^^e"" - + Y^e-' ( - 8^3^ + II7F5 + 27Y,) (3.17) 

- 5 2^/4 e-^'-/' X3 + ^e-^^'- (sV^Yi - 51 y5 + 3 Fe)) 

+ -^2^/^e-^"/^ f - 1188^3X1 + 243v^X3 + 2320X6) + 0(6"^'') , 
585 \ 
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'-'-'4 ^9r/2 I ^< -'4 ^5r/2 



16 23/4m 
3Fi 



+ 



2 
24 
13 



2r, 



26 23/4 

2731^5 + 2A/y; 



(Fi -6V3F6)e- 



2r 



27 X 
10 21/4 

324 
325 

11957859009X3 
28155400 21/4 



3 g3r/2 _^ 350271 Y^^ ^^12 



182872 23/4 m2 



(3.18) 



23/4^^g-r/2^ 



484605^4^^ ^_3,/2 



298792 23/4 



e-5-/2 + 



7978373883 Y^^ 



-7r/2 



21130570240 23/4 



-4r\ 



From these expansions we can extract the UV behavior of the fields which is important to 
understand the holographic physics. For this purpose we have to relate our radial variable r to 
the standard AdS coordinate pAds as 



pAds ~ e' 



3r/2 



(3.19) 



A discussion of the holographic behavior can be found in [IB], where it was shown that the 
integration constants X^ and Ya are paired into normalizable and non-normalizable mode. In 
order to be self-contained we tabulate in Table [l] (which is adapted from [I3]), the leading terms 
coming from each modes. Note that since we obtain the asymptotic behavior from an analytic 
solution, we can relate the integration constants of |il3j to the IR singular behavior of the same 
modes. In particular, one can explicitly check if an IR regularity condition on one integration 
constant is compatible with the absence of the respective non-normalizable mode in the UV. We 
will come back on this point in the next section. In the following table A is the dimension of 
the local operator O holographically associated to the two supergravity modes whose asymptotic 
is (dual to the vacuum expectation value of O) and p^^^ (dual to a deformation of the 
action SS = J d'^xO). Also, the combination which appears at dimension A = 7/3 is the linear 



combination of Yi, Is and Ig which appears in the corresponding terms in (3.16) and (3.17). 



dim A 


non-norm/norm 


integration constants 


6 


PAds/ PAdS 




5 


PAdsl PAdS 


^5 ~ ^e/ -^5 ~ X& 


4 


P AdS 1 PAdS 




3 


P\dsl PAdS 


Yi + F3/X2 


7 
3 


-2/3 / -7/3 
PAdS 1 PAdS 


Y^ + Y^ + YJX^ + X^ 


5 
3 


-4/3 / -5/3 
PAdS 1 PAdS 


l^i/Xi 



Table 1: The UV behavior of all fourteen modes for the S'0(5)-symmetric deformation around 
the CGLP solution, extracted from the asymptotic of our analytic solution. 
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4 Charges and M2— branes 



The space of solutions we solved for in the previous sections should contain the linearized per- 
turbation of the warped Stenzel space due to the presence of a stack of smeared anti-M2 branes 
placed at the tip of the geometry. This configuration was studied in the probe approximation 
in [inj and corresponds in the dual gauge theory to a metastable supersymmetry breaking state. 
In order to identify the backreacted solution, we need to impose the correct boundary conditions 
associated to the presence of the anti-branes at the tip. In this section we start by discussing 
the standard notions of charge in the Stenzel background (see for example [201 Ell 122]) and as a 
warmup we identify the BPS perturbation of the CGLP solution ascribed to the presence of M2 
branes. The anti-M2 brane perturbation will be discussed in the next section. 



In the Stenzel background we can define a "running" M2 charge by integrating -kuG^ on a 



7-dimensional section A^^ = K 2 of the transverse cone at a fixed r 



QM2{r) 



(27r/p) 



4 , 



(4.1^ 



Mr 



where Ip is the Planck length in eleven dimensions. We can also integrate G4 over the 4-sphere 
which has a finite size at the tip and define the quantity 



q{r) 



{2-Klp 



Ga . 



(4.2) 



Si 



For the parametrization (2.1)-(2.2) and for the CGLP background we find from (2.33) (see also 
Appendix IbI) 



6 2" m2 Voly^ 2 



A/2 1 



34 (27r/p; 



6 



^(/.o(/o-2M-^) 



g°(r) 



^^o(r) 



(27r/p)= 



(4.3) 
(4.4) 



where Volv^^ = 277r''/128 ^23]. Substituting the zeroth-order solution (2.15) we find 



Qlnir) 



m? tanh^r 
108 7r2/6 



g°(r) 



m sech r 
3v^7r/3 



(4.5) 



which is the known result for the CGLP solution [TO] . 

We now want to calculate the first-order corrections to these charges from the first-order 
perturbation of the Stenzel geometry. The simpler case is the BPS one, where a stack of M2- 
branes smeared over the is placed at the origin r = 0. The perturbation on the geometry 
should still preserve supersymmetry, so we are forced to set Xa = 0, a = 1,...,6 since the 
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QM2{r) 




r 



Figure 2: The profile of tlie M2 cliarge Qm2 for the BPS perturbation, for different values of the 
parameter Yq. The black dashed line is the zeroth-order solution (Ig = 0). Note that at the 
linearized level the perturbations vanish in the UV. 



"conjugate-momenta" are the modes that parametrize the supersymmetry breaking. Note 
that in this case the solutions for the modes (pa are given by the homogenous solutions of the 
coupled system of ODE's (2.27)-(2.32) and they are easily found by setting = in (2.39). 
The perturbation due to the presence of M2 branes at the tip is found by imposing the following 
conditions on the Ya integration constants: Yi = 1^2 = to cancel IR divergencies in (pi and 
(p2, Y4 = Y^ — Yq = to cancel the divergent terms in the UV expansions of ^4 and ^5, and 
finally Y3 = to fix the freedom of rescaling the three-dimensional coordinates. The first-order 
perturbation to (4.3) is proportional to ho{(p5- 



^e) +/o06 and at the linearized level the running 



M2 charge is given by 

Qm2 = + Qm2 



108 7r2 



tanh^ r + 6^3 Yq (l - tanh^ \ 



The profile of the charge is shown in Figure [2] for different values of the constant Y^. 
asymptotic behavior is the following 



(4.6) 
The 



Q 



IR 

M2 



' Y, + 0{r') 



6^3 



Q 



71^ 



= ^ + C^(e~'0 



The integral of G4 over the four-sphere is given by the behavior of the mode h 



g = gO + gi = — ^ ( 1 - 3^3 Yq ) sech r , 

.3\/.37r V / 



3y37r 

We see that q vanishes at infinity while in the IR it approaches a constant value 



IR 



3^3 



(l-3v^F6) +(^(r'), g^'' = C(e- 



(4.7) 
^6 and thus 
(4.8) 

(4.9) 
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We will denote M = g°(0) the number of G4 flux units through the non-vanishing S*^ at the tip. 
Note that for the zeroth-order solution we have 

(M)2 



4 



Q 



uv 

M2 ■ 



(4.10) 



At first-order, we expect a term related to the explicit brane charge in the IR; in fact, we easily 
see from (4.9) that our solution satisfies 

(g^^)2 



VM2 VM2 



IR 



4 



(4.11) 



which indeed reduces to (4.10) when Yq = 0, which corresponds to having no regular M2-brane; 
Allowing a nonzero Yq introduces a singularity in the warp factor 04 

7 33/4 n 



Hnr^ 



(4.12) 



which is the expected divergency due to smeared M2 branes on the 5*^ at the tip. 

We could have derived these results without relying on the actual solution for the 0^ modes. In 
fact, the linearized BPS perturbation can be obtained by simply shifting the fluxes as follows |13j^ 

05 = 2c , 

06 = |, 

where c is the number of M2 branes. The M2 charge thus changes in the following way 



(4.13) 



Q 



(1) 

A/2 



2 7r2 



hp, 



0<P6 



' ~A 2 

4:71'' 



6^3 7r2 



while for the warp factor we have, from (2.32) 



from which we get 



A0l = cmVo(r) e-3(^«+"°+^«) 
. 33/4 c 



233^ 



(4.14) 



(4.15) 



(4.16) 



which agrees with (4.12) with the identification Yq = c. From this result we can also extract the 
correct mass/charge normalization between the warp factor divergency and the charge sourced 
by the branes, which will be useful in the next section 



(4.17) 



^Note that equation (4.11) is just the standard relation introduced in |2(J]. The quantities of equation (2.16) 
of that reference are $ = Q^^^, N = Qif^ and J GAG - 

"^By shifting the fluxes f-^f + 2c, h^h 
terms proportional to which are not seen in our linearized deformation space. 



I)- 



^ (see also 

I we can obtain the full nonlinear solution, but this introduce 
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In the next section, we will turn to the case of interest in which we add a stack of anti-M2 



branes at the tip of the transverse Stenzel space. In this case the expressions (4.3), (4.4) evaluated 
at first-order in perturbation theory will depend on all of the Xa, Ya integration constants. 
However, we have to impose appropriate regularity conditions for the IR and UV behavior of 
the modes (pa, and we will see that this fixes all the integration constants in terms of X4, which 



is the one responsible for the force on a probe M2 brane in this background (see section 5.3), 
and Yq. We thus expect the expressions for the charges in the BPS case to be modified by some 
pieces proportional to X4. By requiring the variation in the M2 charge Q^/L t)e commensurate 



to the singularity introduced in the warp factor, equation (4.17), we will derive a relation which 
fixes Yq in terms of X4 and so we will fix all the integration constants in terms of the number of 
anti-M2 branes. 



5 The anti— M2 brane perturbation 

In this section we consider the perturbed solution corresponding to a stack of N anti-M2 branes 
at the tip of the transverse geometry. It was shown in pll] that in the probe approximation, for 
N/M < 0.054 this configuration is metastable and will eventually decay into a supersymmetric 
configuration in which M — 1 — N M2 branes are present at the tip In order to find the 
supergravity dual to the metastable state, we will adopt the following strategy. Firstly, we 
consider the IR behavior and we allow only for divergencies that are directly sourced by the anti- 
branes. Secondly, we demand that the UV non-normalizable modes described in section 3.3 are 
absent, so that the UV asymptotic is the same as for the original CGLP background. As we will 
show, these requirements (together with the mass/charge normalization discussed in the previous 
section) provide enough independent contraints on the deformation space to fix every integration 
constants in terms of a single physical quantity, namely the number of anti-M2 branes present at 
the tip. We then compute the relevant charges for the perturbed solution, as well as the explicit 
expression for the force felt by probe M2 branes in the backreacted anti-M2 background. 

5.1 IR and UV boundary conditions 

We now proceed to impose regularity conditions on the IR behavior of the modes (pa- We demand 
that divergencies are zero except for the singularity in the warp factor 04 which is directly sourced 



The units of G4 flux for the susy state are then M — 2. A way to understand these values is to look 

,uv _ 

susy 



at (4.11). We then see that these are the correct values so that the charge at infinity is conserved: Q,^^ 



- — + M — 1 — iV— = ^ — = , where Qsusy and are the charges for the susy and metastable states. 



18 



by the anti-M2 branes. We first impose the zero energy condition, which amounts to setting 



Xo = 0. 



(5.1) 



From regularity of </)i we derive 
2 



Xi = - 



X, 



8 V3(X5 + Xe) + 15 X4 Ho 



27 

27^3X3 SXg 



(5.2) 



m^X^ (lOOi/o^(-l) + y2 3^/^A;^J 



20 5 20y2 3i/4i^(-l) 80y3ir(-l) 

while from the singular terms in <^2 we derive Xg and Y2 in terms of X3, X4 and Yi 
9v^X3 Fi m2X4 (220i7oi^(-l) + y2 3i%J 



Xfi 



12y2 3V4/s:(-l) 



48v/3ir(-l) 



(5.3) 



V-594.„i/4 v^3V^(9 3V^ + go)yi m^X4 

^'-^^^^ 5K{^) lOi^(-l) 

+ 30K{-l)k2 + 8ifo(153 3^/^ + 2Ho)K{-l] 



We can check that with these conditions the other IR divergencies of the modes (pa are automat- 
ically canceled, except for a logr mode in the IR behavior of the field 03, which is a perturbation 
of the metric, which is proportional to X4. It is not clear why one should not be able to kill 
such divergent behavior. However, after imposing the previous boundary conditions, the solution 
presents an even worse singularity appearing in the field strength F| 



p2 

4 ^45 



(5.4) 



which is quite analogous to the divergence found in the anti-D3 solution [7], with the difference 
that now the action is divergent. This behavior is sub-leading with respect to the energy density 
associated to the divergency in the warp factor, which is of order r~^. Note that this is an IR 
phenomenon insensible to UV boundary conditions; in fact, the integration constant X4 cannot 
be set to zero for the very simple reason that it parametrizes the force felt by a probe M2 
brane |T3] and thus is indicative of the presence of anti-M2 branes at the tip. Despite arguments 
in the literature, there is not a rigorous proof that shows if this singularity is acceptable or not|^ 
Given the difficulties in proving this, we will assume that the singularity is harmless and we 
will try to see if the anti-M2 solution develops problems in the UV; if this is not the case, the 



^For the anti-D3 case, it was argued that the singularity is an artifact of perturbation theory and will disappear 
in the full backreacted solution (see [21 [25]). However, other works pointed out problems in the full solution for 
antibranes in ISD flux backgrounds [261 El] ■ 
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solution we find describes the holographic dual of the conjectured metastable state in the field 
theory, but clearly a more detailed study of the IR singularity is required to decide whether this 
supergravity solution can be trusted or not. 

We now proceed by imposing boundary conditions in the UV, where we demand that non- 
normalizable modes in the UV expansions for the modes (pa are absent. The first condition is 
from the e^^^^ term in 02? from which we get 

X3 = 0. (5.5) 

From the divergent term in 05 we get 

n = , (5.6) 

and finally from the term e~^^ in 0i we get 

ri = 0. (5.7) 

Note that we should allow an e"'' term in the fluxes, which is dual to the dimension A = 7/3 
operator, since it is the charge mode sourced by the branes. We thus see that we fixed the ten 
physically relevant integration constants in terms of X4 and Y^, which are related respectively to 
the force on a probe M2 brane and to the number N of anti-M2 branes placed at the tip 

5.2 Charges and anti— M2 branes 



In order to relate X4 and Yq we look at the M2 charge (4.1). Once all the boundary conditions 
are imposed, we get that 

where the coefficient a is the following combination of the numerical constants which enters in 
the expansions for the modes (pa 



a ^ ^°'°^^'^^ + ^^^°' - + + f ^°>**- ^ 0.900178 . (5.9) 

60^3 ?,m^2K{-l] 



We impose that this variation gives the correct singularity in the infrared expansion for the warp 
factor, which is found to be 

33/ny.-^..-A,) ^ (5.10) 

with 

/3 = a + ~ 0.958828 . (5.11) 
v3 
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From the mass/charge normalization (4.17) we thus get the following condition 



(5.12) 



which results in 



Ye = ^ia + P)m^X4= (a + X^. 



2v/3/ 



(5.13) 



If we now plug this relation back into the expression for the charge (5.8), we find the following 
relation 

Hi X4 



Q 



IR 
A/2 



36 7r2 /6 



(5.14) 



We note that this result does not depend on the UV boundary conditions. Indeed, although it is 
not clear from our derivation, it is easy to show that if we only impose IR boundary conditions 



the terms proportional to X3, Y5 and Yi that appear in (5.8) and (5.10) cancel in (5.13) 



Since Qif2 is related to the number of anti-M2 branes placed at the tip, from (5.14) we 



determine X4 as a function of and thus we fix all the integration constants in terms of this 
parameter. 



With these results, we can explicitly compute the charges associated to the anti-M2 brane 
perturbation (in Figure M we show the profile of the first-order perturbation to the Maxwell 



charge (5^/2)- particular, the M2 charge Qm2 evaluated at a holographic screen at infinity 
should be the same for the metastable and the supersymmetric state. This condition ensures that 
the metastable state is a state in the same theory which is dual to the supersymmetric vacuum. 



Unfortunately, we see from (4.7) that the perturbation to the M2 charge vanishes at infinity at 
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the linearized level, and so in our backreacted solution the value of fixed. We expect 

shifts of this quantity to appear only at second-order in perturbation theory. While we cannot 
directly check whether the value of the charge at infinity is conserved, we can look in the IR and 
check whether the charges are perturbed in a consistent way. From the probe computation (see 



footnote^, we expect relation (4.11) to be satisfied. For a supersymmetric domain wall, this 



easily follows from the equation of motion (2.4) and the self-duality of the fiux G4, and indeed 



we found that the BPS perturbation considered in section |4] is consistent with this constraint at 
the linearized level. For the non-supersymmetric case, one should be more careful. It is useful 
to write the first-order perturbation to the Maxwell charge in the IR in the following way 



Q 



IR 

M2 



27r2 



1 



6^3 



05(0) -406(0) +—=06(0) 



3^3 



(5.15) 



from which we derive, at the linearized level 
{qiR? 1 



6^3 



+ 06(0) 



Q^i2 Q 



IR I 

M2 "T 



12^3 



7r^ 



05(O)-406(O) +0{X'). (5.16) 



After imposing the anti-M2 IR boundary conditions, we find that the term in the brakets in the 
right hand side of the last equation is not zero 



5(O)-406(O) = /fo3-^/'m2X4. 



(5.17) 



Indeed, this is the term which gives rise to the singularity in the filed strength F| that we 
discussed in section 15.11 This result is consistent with the fact that at the linearized level the 



self-duality of the four-form fiux is spoiled, and we do not expect relation (4.11 ) to be satisfied 



for the anti-M2 solution. As we discussed in the previous subsection, it is possible that this 
result is an artifact of the perturbation theory. While we cannot address this issue within our 
first-order technology, we believe that further investigation is needed on this problem. 



5.3 The force on a probe brane 

With the results obtained in the previous subsections, we are able to compute explicitly the 
coefficient of the force exerted on a probe M2-brane in the anti-M2 backreacted background, 
whose functional form has been derived in flSl [T8] : 

ISX.cschV 



M2 



(2 + cosh2r)3/4 ' 



(5.1^ 



Inserting the expression for X4 that we derive from (5.14) we obtain 



M2 



648 TT^ /^XcschV 
(2 + cosh 2r) 3/4 



(5.19) 
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This result has to be compared to the one given by the probe anti-brane potential a la KKLMMT 
[28], which is given in [18] and reviewed in Appendix pj The result of this computation is given 



in (B.IO ). Once we substitute c?2 we see that the two expressions exactly agree. This is a nontrivial 
check that our IR boundary conditions are the correct ones to describe anti-M2 branes in the 
Stenzel geometry. 



5.4 Asymptotic of the anti— M2 solution 

We now collect the results we obtained for the twelve (Xq, Ya) integration constants and which 



determine the anti-M2 solution in terms of the constant X4, which is fixed in terms of by (|5.14|) 
For the Xa integration constants we have 



36 7r2/6 _ 
N 



(5.20) 



-2Hom^Xi, X2 = 0, X3 







(5.21) 



-2 



m2X4 
12^3 



12^3 



For the Ya integration constants we have 



Fi = 
Y2-- 

>3 = 0, 



(5.22) 



-2 H 



4 , 



v^/7o(633i/^ + 52/7o) -180fc^^ +2(27 + 3^/4 /Jo)v^A;^,r Q 



m2X4 
180 



The IR and UV behavior of the backreacted anti-M2 solution, up to the desired order, can 
be read off from the analytic solution presented in Appendix |X] after imposing the boundary 



conditions (|5.21|), (|5.22|). For the reader's convenience, we show here the first few terms of the 

(5.23) 



ultraviolet behavior of the solution. 

~ 4m2X4 ' 
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IGm^X^ 
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29 2^/^HoK{-l) + 6^/^ k^^ 
1102^/^HoK{-l) + 6^/'^ k^^ 
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-9r/2 



23 



+ 



3i^(-l) 
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Qm^Xi 
m^X4 
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-9r/2 
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4244 23/4 1466 2^/4 A;t 



1755V^ 3315 31/4 ii'(-i) 



,-lW2 + C»(e-7r) 
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+ 



65 

2192y3m2X4 
14365 



V3Ho{63 3^/^ + 52Ho)-180k^^+2{27 + 3^/^Ho)V^k^J 



-2 



-2r 



^3 Ho{Q3 3^/^ + 52 Ho) - 180 k^^ + 2 (27 + 3^/^ Ho)V^ k^^ T (j 



-4r 



26172=^/^ iJo^(-l) + 80 6^/^ k^^ 



6r\ 



(5.28) 



6 Conclusions 

In this paper we constructed the analytic solution for the twelve-dimensional space of linearized 
non-supersymmetric deformations of the warped Stenzel space, consistent with the S0{5) sym- 
metries of the supersymmetric background. Our solution provides an interpolation between the 
IR and UV behaviors previously constructed in [13j and it should contain interesting informations 
about the dual (2-|-l)-dimensional gauge theory. In particular, we were interested in finding the 
supergravity solution dual to metastable states, which were conjectured in [10] to be described 
by a stack of anti-M2 branes placed at the tip of the transverse geometry. We were able to 
identify this solution by imposing suitable boundary conditions on the set of twelve integration 
constants (X^, Ya) that parametrize the full deformation space, and indeed we showed that this 
solution is unique and it depends only on the number N of anti-M2 branes placed at the tip. We 
then used this solution to compute the force exerted on a probe M2 brane placed in the anti-M2 
backreacted supergravity background and we showed that it exactly agrees with the calculation a 
la KKLMMT [25] in which one considers the anti-M2 brane as probing the backreacted geometry 
of M2 branes on the Stenzel backgrounc ^ 

The linearized supergravity solution displays however an IR singularity in the four-form flux, 
which leads to a divergent action, whose nature is still poorly understood. Our analysis shows 
that this is the only drawback of the supergravity solution, which otherwise has the desired 
features to describe the metastable state in the dual gauge theory. It is thus of great importance 
to establish the nature of this singularity. However, proving if these singularities are acceptable 
or not is unfortunately beyond the reach of our first-order technology and clearly more work is 
needed to address this problem. 
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A Analytic solutions 

Here we show the analytic solutions for the modes which can be obtained by explicitly per- 
forming the integrations that appear in (2.35) 

I4 = m^X,H{y) , (A.l) 

^5 = -2V2{y' - 3)-\y' - l^Uy) - T'l\y' - 3)-\y' - lf'L,{y) , 

Is = ^V2iy' - 3r\y' - ly'^'iSy' - ^)L,{y) - V2{y^ - 3)-\y' - l){y' - lf''L,{y) , 

6 = \x,{y' - 1) + \x,{y^ - V + 3)^/^ + ^(y^ - 3f{y' - l)LM + ^{y' - ly^LM 



4 



{y' - 4)L,{y) + m'X, (^(|/^ - l)H{y) - 2V2y-%y'' - ly'^^y' - 3) 



where 



3^3 

L,{y)=X, + m^X,G{y) 



(A.2) 



L,{y) = Xe - m'X,(G{y) + ^H{y) 
Lsiy) = Xs + m X^i^ 



22G{y) UH{y)^ 



.27y3(y4_3)(^4_ 1)3/2 27v^ 

We recall that the variable y is defined as 

y={2 + cosh(2r))^/^ 



27 J 



(A.3) 



(A.4) 



and the expression for the warp factor H{y) and the Green's function G{y) are given in (2.16) 
and ( |237l ). 

We now show the expanded form of the solution for the modes (pa in terms of the variable 
y, obtained by replacing the analytic solutions for the modes C,a in (2.39). We impose the zero 
energy condition, so we put X2 = 0. 



1 



(3-4|/4 + J/8)l/2 



27 



V2{y^ - 3)1/2 ^2 



^ (y4-l)V3(y4-3)V2 
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^(^4 _ 3)1/2 ^33/4 

1 2v [ r 3v^(-19 - 26w^ + l^u^)H{ 
+ 2^^^; - 3)3/2 



- 1)2 



- (y' V + 3) V2 + ( ^^^^^"(^') I ^) - 32V3F(arcsin(,^)|i) 



(A.5) 



- 301 4 

02 = ^ + 



^4 ^4(^4 _ 3)2 



+ 9^^^V^ + 9^3yv^ (825 - 63^ + 343,^ - 49,-) 



140v^ 



+ ^'^^4'!.^i)3/2 + ^^^(63Xl + 432X3 - 7v^(5X5 + llX6))F(arcsm(y-^)| - l) 

. 7 o 2^ f r (-171 - 342x^4 + 936w« - 546^^ + %lv}^)H{u) , 

- 4V6X6I/ _ 1 + 3771^X4 ( / ^ ' ^ ' 



12v^(li4 - 1)1/2 

y {u^ - 3)2(99 - M2u^ + 176ii« - 154x^12 + 17u^^)G{u) 



6a/6(m4 - 1)5/2 



, 2t/2(138-119y^ + 14t/S) , 16 2 
+ ^(I^^i^ + yarccoth. 



-du 



(A.6) 



03 = 1^3- 2) 01--02 + 



3 ~ , 21X^y^^f^ , 27X31/^^^ , v^Xgi/lSy^ - 7) 



+ 



+ 



8"^ ■ 8^2(1/4-3) ' 4v^ ' (y^- 1)^/2 

+ ^ (^135X1 + 432X3 - 8^3(5X5 + 12X5)) F arcsin(y-^) | - 



+ ^ (7v^Xi- 4(^5 + X6))g(|/) + 3771^X4 j 



+ 



J' (-33 + 1257i^ - 797/« + llu^^)G{u) ^ , St/^ V3_„_,_„,_ / y 



2v^(774 - 1)5/2 



-du + 



;-51 + 9m^ + 39^8 - 13m-)/J(m; 

4V^(k4 _ 3)3/2(3 _ 4^4 + ^8)1/2 

1 + 



arctanh ^ + log . 
4-4y4 4 VV3/ Vl-y 



(A.7) 



For the flux perturbation we have 



,V2V. V. -My'-^w-^y'" 



\ ( \ 

A5 



(A.8) 



27 



where 



A5 = n + 



6v^ 



+ 



5v^(l/4 - 1)3/2 



10v^(l/4 - 1)3/2 



(A.9) 



1072X5^ 



9(y4-3)(y4_ 1)3/2 



+ 



272X67/(272 - 27%^ + 60?/^ 
45(2/4 -3)(?/4_ 1)3/2 



108 



27Xi+44y3(X5 + X6))G'(y) 



1 / 120fX + X ) \ 

+ 180 ( " ^^^^1 " 297 V3X3 + 40X5 + lOeXg - SOXgy^ + — j i7(^) 



+ 771^X4 



Hiu) 



76 + 85u* - 78m« + 13^12 2 3 ^ , ^^^^ ^ 8y^H{y) 



2v^(m4- 3)3/2(^8 _ 4^4 + 3)1/2 3 



77^ + ;TrG(y) + ^H{y) + 



+ 



^ du H ^ arctann 

3V2{u^ - 1)5/2 ^4 _ 1 3 



y 



3 

2 



V3{y' - 3)= 



log 



1 + y^ 



, ^ y\-S + y'f ~ ^ 3v/3Xi|/(7-5y4) 3v/3X3y(72 - 45/ - Sy^) 2v^X52/(/ - 3) 



6v^(-l + 7/4)2^^ 572(^/4-1)3/2 



10v^(y4 - 1)3/2 



3(y4 - 1)7/2 



^ 2v/2X6y(20y4 - 33) ^ ^^^^ / _ 7x/3Xi _ 9v/3X3 ^ X6(ll - 5^^) ^ X5(l + 6^^ - 32/8) 



15(2/4 - 1)3/2 



40 



10 



30 



6(2/4 - 1)2 



+ m^X4 



Hiu) 



2/4(-19- 262/4 + 132/8) 2y3(/ _ 3)(2/8 + 3)^(2/) , 2/^^(2/) 



2v^(2/4 - 1)5/2 



+ 



3(2/4 - 1)3 



+ 



M4(y4_ 3^2^^^ _ 38^4 + ^1^8)^(^) ^ 22/='(-3- 42/^ + 32/*) ^ 1 



3^2(^4 - 1)9/2 



3^3(2/4 - 1) 



x/3 



-1-2/2 



B Brane/antibrane potential 

In this Appendix we review the calculation of the force on a probe antibrane in the Stenzel 
background which has been performed in [18]. We consider a stack of M2 branes at a position 
r = ro in the transverse geometry and we want to compute the force exerted on a probe anti-M2 
brane placed at the tip r = 0, due to the backreaction of the M2 branes. To compute the full 
backreacted geometry we only need to add a harmonic function 6H{r) to the background warp 
factor Ho^r) [30]. This function is given by the Green's function on the warped Stenzel space [29] 
and since we are considering smeared branes, we only need to solve the radially symmetric Laplace 
equation. The Laplacian is given by 

1 d f.r7^j^d5H\ 



^SH = ^TTT VG(/""— , (B.l) 
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where we are labeling the eleven dimensional coordinates by x/ (/, m = 0, . . . , 10) and G = det gim- 

(B.2) 



From (2.1) we easily get 



G 



gvr ^ g-2-27 



and so imposing ASH = we get the following equation for 6H'{r) 



const 



(B.3) 
(B.4) 

(2 + cosh2u)3/4^' ' ^^'^^ 

and we should set 6H = Hi for r < tq and 6H = H2 for r > tq. The constant di is then fixed 
from the matching condition di = H2{ro) and the constant d2 is related to the number of M2 



So the two solutions of the Laplace equation are, using (2.15) 

Hi = di 

csch 



Ho = d, 



dr . 



branes from (4.1). In fact, we have 

N - I - ^''"^'^"^^-'^ H' e^(-+/^o) 



(B.6) 



where 5M. is a small shell around tq. If we use H2 given by (B.5), the above equation thus fixes 
d2 in terms of N 

(B.7) 



3247r2/-6Ar 
do = 



We now compute the force exerted on the probe anti-M2 brane by looking at the variation 
in the potential when we move the stack of M2 branes away from r = tq. For anti-M2 branes, 
Vdbi = Vwz and since we have 

1 



V, 



DBI 



-32 



im? H ' 

the potential is just proportional to 2H~^ . Expanding this we get 



V 



m? H m?HQ 



6H 
Ho 



(B.8) 



(B.9) 



and so we easily get the force: 



2 dV 



' M2 



m? dro 
1 



r=0 



2 d2 csch 



m2i7o2(2 + cosh2r)3/4 ■ ^^'^^^ 

This result agree with the computation of the force exerted on a probe M2 brane due to the 
backreaction of a stack of anti-M2 branes (5.19). 
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